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Abstract. This paper deals with classical solutions of the modified chiral model on 
R 2+1 . Such solutions are shown to correspond to products of various factor which we call 
time-dependent unitons. Then the problem of solving the system of second-order partial 
differential equations for the chiral field is reduced to solving a sequence of systems of 
first-order partial differential equations for the unitons. 
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I. INTRODUCTION 



As is well known |2|], some physical phenomena (relativity, particle physics, etc) 
lead, in their mathematical description to nonlinear partial differential equations. Some 
of these equations are elliptic (when we deal with time independent problems), others 
are hyperbolic or parabolic. However, with few exceptions, nonlinear partial differential 
equations are hard to solve (especially when they involve more than one spatial vari- 
able). Although in some cases static solution can be found the generalization to the time 
dependent problems often presents a formidable task. 

A simple model which exemplifies many of the points mentioned above is the so-called 
relativistic chiral model in (2+1) dimensions. This model involves a chiral field J(x,y,t) 
which takes its values in SU(N). The equation of motion for this field is given by 

d^J^d^J) = 0, JeSU(N), (1) 

and is, in fact, the Euler-Lagrange equation corresponding to the Lagrangian density 

C = ^tT[(J-%J)(J-^J)}. (2) 

Here the indices }Ji,v,a range over the values 0,1,2, with x^ = (t,x,y) and <9 M denotes 
partial differentiation with respect to x^. Thus J describes maps of R 2+1 into SU(N) 
and, when we restrict our attention to static fields, we have maps from R 2 into SU(N). 
If we restrict our attention further and consider only static fields whose energy is finite 
(see below) then the problem reduces to finding harmonic maps of S 2 — > SU(N). 

Such harmonic maps were studied by Uhlenbeck || who, in a seminal paper, showed 
how to construct these maps in an explicit form. Her construction involved the intro- 
duction of the concept of a uniton. She then showed how starting from a given solution 
one can find a new one by "the addition of a uniton" . This observation allowed her to 
demonstrate that all solutions can be obtained from the trivial ones (ie J=const) by the 
addition of a certain number (in fact < (N — 1)) of unitons. 

This beautiful result relied heavily on the integrability of the chiral model in (2+0) 
dimensions (for more details see PJ). When we consider the same model in (2+1) di- 
mensions, ie when we want to solve ([]]) we immediately run into a problem as it is not 
integrable and so Uhlenbeck's method does not apply. 



Some time ago Ward || presented a model in (2+1) dimensions which is an integrable 
generalization of the static chiral model in (2+0) dimensions. The model involves an 
SU (N) valued chiral field whose equation of motion now takes the form 

(Tf» + e> u ' a V a )d lt (J- 1 d v J) = 0. (3) 

Here the tensor i] tiu is the inverse Minkowski metric, given by if" = diag(— 1, 1,1); s^ ua is 
the fully antisymmetric tensor with e 012 = 1; and V a is the constant unit vector pointing 
in the x-direction, V a = (0, 1, 0). 

Notice that although (^j is not relativistically covariant (as V a is nonzero) the solutions 
of Uhlenbeck are static solutions of (|3|). As the model is integrable many solutions of (|3|) 
can be easily found. The first solutions of this model |5| exhibited many properties 
seen in other integrable models (like wave-like behaviour with a phase shift, soliton-like 
structures passing through each other, etc). Recently, Ward |J, Ioannidou and even 
more recently Anand |§ have found further interesting solutions of (§). These solutions 
described the evolution of soliton-like extended structures whose properties were quite 
similar to properties of solitonic objects in the 0(3) sigma model (ie they exhibit the 90° 
scattering). 

The approach of Ward and Ioannidou involved treating fields of the model as the 
limiting cases of the Riemann problem with zeros. In this paper we return to Uhlenbeck's 
construction of static solutions and show that her construction can be generalized to give, 
in general nonstatic, solutions of the modified chiral model of Ward (ie of (0)). Our 
approach reproduces the solutions reported in and and presents an unusual but 
interesting application of Uhlenbeck's approach. 

II. THE MODIFIED MODEL 

As shown by Ward |JJ the modified chiral model ([|) has the same conserved energy 
density as the original chiral model (0) and this energy is given by 

S = -\ t r[ { r^f + (J-^f + { J- l Jy)% (4) 

To specify the problem of finding solutions of (^ completely we have to state the 
boundary conditions satisfied by J. Thus in order to ensure the finiteness of energy of 



the solutions, we require 

J=J c + J e (9)r- 1 + 0(r- 2 ), r -> oo, (5) 
where z = re w (with z = x + iy), J c is a constant matrix, and Jg is independent of t. 

Equation ([3]) is completely integrable, in the sense that can be written as the compat- 
ibility condition for a pair of linear equations involving a spectral parameter A. Define 

A z = \j~\d z - l -id t )J, A- z = \j~\d- z + Ud t )J, (6) 

where A z , A z are N x N matrices depending on x, y, t but not on A. 

Then the pair of linear operators, 

0, + (1 - X)A- Z , 0, + (1 - A- 1 )^, (7) 

with Z = d z + |zA _1 <9 t , 0g = d z — \ i\d t) commute for all AeC* = C — {0}, as a 
consequence of @. Therefore one can find a solution of £a of the equations 

AE X = (1 - A- 1 )^^, tf 2 -£ A = (1 - A)^ 2 - (8) 



Theorem 1 If J : R 2+1 — > SU(N) is harmonic, then there exist a GL(N,C)-valued 
function E\ of (§) which is analytic, holomorphic in AgC* and unitary for |A| = 1 or 
A = A -1 suc/j t/iat ELi = J and 

Ex = (E^y, (9) 

where > denotes complex conjugate transpose. Conversely, if E\ is holomorphic and ana- 
lytic; and A z , A z defined by (|6|) are independent of \, then J = E_i is harmonic. 

This follows from || and E\ is called an extended solution corresponding to J. System 
(§) is overdetermined and in order for a solution E x to exist, A z and A z have to satisfy 
integrability conditions, 

(d z + ~id t )A z + (d z -~id t )A z = 0, 

(d z -~id t )A- z -(d z + ~id t )A z -2[A- z ,A z ] = 0. (10) 
Notice that, the first equation of fllOD is equivalent to (||). 



III. UHLENBECK'S CONSTRUCTION 



It has been known for some time that all static solutions of grassmannian models are 
also solutions of the static chiral model |1; although, not much has been known about 
other solutions. However, Uhlenbeck [|J proved that any harmonic map from the Riemann 
sphere S 2 to the unitarity group U (N) can be factorized into a product of a finite number 
of factors (so-called unitons) involving maps into Grassmannians Gk{C N ) each of which 
satisfies a system of first-order partial differential equations. Thus the problem of solving 
a system of second-order partial differential equations for a static harmonic map has been 
reduced to having to solve a sequence of systems of first-order partial differential equations 
for the unitons. This observation can be generalized to the non-static case (0). 

We will seek solutions of (||]) of the form 

n 

Ex= l[{l + {X-l)R k ), (11) 

fc=0 

for some n, where Rk are Hermitian projectors of the form Rk = (q\, <g> qk)/\qk\ 2 (Ro is a 
constant matrix, ie 0-uniton), which satisfy some first order differential equations; and q^ 
are iV-dimensional vectors which, in general, depend on z,z,t and \qk\ 2 = qu • %■ Note 
that a constant E\ is a 0-uniton. Then, following Uhlenbeck, we consider 

E x = E x {\ + {\-l)R n ), (12) 

and find 



Theorem 2 Let E\ : C* xR 2+1 — > GL(N, C) be an extended harmonic map E\ : R 2+1 — > 
SU(N) for |A| = 1. Then Ex = E\(l + (A — l)R n ) is an extended map for a Hermitian 
projection (R n ) if and only if 

(1 - R n )(d t R n - 2iAzR n ) = 0, 
R n {d- z R n - A- Z {1 - R n )) = 0, (13) 

where A s = — (A z y given by (||). 

From the definition of the extended solution, A 2 = (1 — A) -1 ^ 1 $zE\ and so 

A- z = (1 + (A" 1 - l)R n ) U- Z {1 + (A - l)R n ) - (d- z - X - iXd t )R^ , (14) 



using the fact that E x l = (1 + (A -1 - l)^)/^ 1 and A z = (1 - X^Exd.E^ 1 . But, from 
Theorem 1 we know that E\ is an extended solution if the A and A -1 terms vanish, which 
leads to the system fllTf). 

Next we observe that 



Theorem 3 Suppose that E\ is an extended solution and so is E\ = E\(l + (A — l)R n ). 
\El x d- z E x and A- z = \E X 



Then if A z = \E X 1 fi^Ex and A z = hE x 1 $ z Ex, they are related by 



A- z = A- z -d- z R n + Ud t R n . (15) 



This follows automatically from the terms in ([LJ) which do not contain A. In general, by 
induction, one can show that 

1 n 

A- zn = {-d- z + ~id t )Y,Rk, (16) 
and so, that the corresponding Rk satisfy the following equations 

n— 1 

{l-R n )(d t R n -2iJ2(-dz + ^id t )R k R n ) = 0, 

fc=0 1 
n— 1 i 

R n (d- z R n -Y l {-9- z + 1 ,idt)R k {l-Rn)) = 0. (17) 



,<■: I) ^ 



Observe that for n = 1, system (|TTD gives d t Ri R\ — and R± d z R\ = 0; ie the equations 
for the holonomic static solutions (instantons) of the grassmannian models. This is the 

case as 

odd- f dtqj g gi _ t gl g gi n t-x o\ 

d t Ri Ri = qi q[ — p-rr d t q{ q x = 0, (18) 

ml ml 

which is satisfied, if and only if, q\ is independent of t. Using the same argument for 
R\d z R\ = we deduce that q% = qi(z). For n ^ 1, we have more general and, in general, 
t-dependent solutions. 

Now, let us concentrate on (|j) and discuss the construction of its solutions. If we set, 
E-i = J and ELi = Jo in (|T2|), then the chiral field J (due to Theorem 2) is of the form 

J=J (l-2R n ), (19) 

and satisfies (§) if J does and if R n satisfies ([13]). This describes a time-dependent 
multi-soliton solution of @ since it satisfies fllCf ), ie 

(d z + l -d t )A- z + (d- z - l -8 t )A z = (d z + % -d t ) (-(1 - 2R n )(d- z + l -d t )R n + (1 - 2R n )A z {\ - 2R r , 



+ (9-z ~ ~3t) ((1 - 2R n )(d z - l -d t )R n + (1 - 2^)4,(1 - 2R ri 



(d z + l -d t ) [A- z - {d- z - l -d t )R^ 
+ (dz ~ \d t ) [A z + (d z + % -d t )R^j 



using equations (|), (|13|) and (|19|) . 

So, to construct time-dependent solutions of the SU(N) model (|J), one can start 
from a constant solution (zero uniton) and add to it one uniton. This solution will be 
nonconstant - but it is static. Then, to this one uniton, one can add a second uniton, 
then a third one, and so on. As we are going to see, the number of soliton-like structures 
in a n-uniton configuration is not related to n {i.e. the procedure of the "addition of a 
uniton" can both increase or decrease the number of unitons). In the static chiral model 
(H) Uhlenbeck's theorem shows that for U(N) the largest uniton number is less than N. 
In our case, the t-dependence of the system (|17|) invalidates this argument; the number 
of unitons can be arbitrary. In the static SU(2) case considered by Uhlenbeck the only 
solutions are those described by constant matrices (0-unitons) and factors constructed 
from holomorphic functions (1-unitons). In the modified chiral model the 0- and 1-uniton 
solutions are the same but then, in addition, we have further solutions corresponding to 
two and more unitons. These additional solutions are nonstatic. We do not know, at this 
stage, whether there is any bound on the uniton number so that all solutions correspond 
to field configurations of up this number. 

IV. THE SU(2) AND SU(3) CASES. 

In this section we shall construct and discuss, as an example, a 2-uniton solution of the 
SU (3) model @. First of all, observe that all fields of the SU (2) model can be embedded 
into the SU(3) model so that the SU(2) solutions are automatically also solutions of the 
SU (3) model. In what follows we will recover some of them as special cases of the derived 
expressions. 

Let us observe that the following expression describes a non-static 2-uniton solution: 

^-Hi^Hi^)' <2i) 



where qi for i — 1 , 2 are the three-dimensional vectors given by 
Qi = (!> /> 

g 2 = (l + \f\ 2 + \g\ 2 )(l,f,g)-2i(tf' + h)(f,-(l + \g\ 2 ),gf) 

-2i(tg' + h 1 )(g 1 gf,-(l + \f\ 2 )), (22) 

and K is a constant matrix. Here / and g are rational meromorphic functions of z. The 
field J takes values in SU(3), is smooth everywhere in R 2+1 (as the two vectors are 
nowhere zero). It satisfies the boundary condition (Q) and the equation of motion ([|). 

Now let indicate how the solution (^l|) was constructed and how further solutions can 
be obtained. One way of proceeding is to take the qi vectors of the form = (l,Wj,t>j). 
Then, for the 1-uniton solution corresponding to (|17|), the qi vector q± = (l,Ui,Vi) (as 
we have shown) is given in terms of the complex variable z. Thus we write q\ = (1, f,g) 
where / and g are functions of z only. Then the system fllTD for n = 2 reduces to the 



following first order differential equations for the function u 2 and v 2 : 

G G 

d- z u 2 = , I j| 2 . I , 2 A -u 2 + f), d- z v 2 = , „| 2 , , 2 A ~v 2 + f), (23) 
(i + l/r + M) (1 + l/l+M) 

where G = /'(—u 2 + /) + ^(-^ + g) + («20 - v 2 f)(fg' - f'g) and 

to = 2? I f 12 i I l2 A Mff +gg)+g{f g-gf) + f), 
(i + l/l + l#l) 

to = 2 * [i * + Mff + - /(/g - 9'f) + </), (24) 

and ' denotes the derivative with respect to 2. A few lines of algebra allow us to show 
that q 2 of ( P2"|) is the general solution of the equations of the above system. Notice that 
when Vi = (ie g = 0), ( |2T| ) corresponds to the field J of the SU(2) model @. 

Calculating the energy density for ( pip we find, 

£ = £ Q + 4tr((d z - t -d t )R 2 (d z + i -d t )R^ -Atx[{d- z - l -d t )R 2 d z R 2 R^ 

-4tr [{d z + ^ifc R 2 d- z R 2 ^j , (25) 

where So is the energy density of the 1-uniton (ie static) field. As is well known |J £q is 
a total derivative and so, probably, also is £ but, so far, we have not managed to prove 
this. However, we have looked at some special cases and have found that the total energy 
E, obtained by integrating £ , is quantized in units of 8ir, ie is given by E = 8Nit. N 
appears to be related to the number of topological structures. 



Let us discuss some special cases: 



(a) The embedding of SU(2) solitons, where we set 

g = f, h = h x , (26) 

Such soliton solutions are localized along the direction of motion; they are not, however, 
of constant size. Their height (maximum of £), is time dependent. In fact, for (deg/ > 
deg h) no scattering occurs, ie the solitons are located at the centre-of-mass forming a 
totally symmetric ring configuration. Otherwise, there are (deg / — 1) static solitons at 
the centre-of-mass of the system accompanied by N = deg h — deg / + 1 solitons (located 
at t + z N = 0; since J departs from its asymptotic value J c when tf' + h = 0) accelerating 
towards the solitons in the centre, scattering at an angle of n/N, and then decelerating 
as they separate. This phenomenon has been observed, firstly, in the SU(2) model (|3]) 
(cf. H [7|) and is a feature of systems (with nontrivial topology) admitting topological 
solitons. 

As an example, let us take f(z) = z and h(z) = z 2 . The corresponding field represents 
a configuration of two solitons that undergo a 90° scattering. The energy density of this 
solution is 

£ = 32 1 + 12r 2 + 12r 4 + 32r 2 t 2 + 8t 2 - 16(x 2 - y 2 )t 
[l + 12r 4 + 4r 2 + 16(x 2 -?/ 2 )t + l + 8t 2 ] 2 ' 
and is symmetric under the interchange t — > —t, x —>■ y, y —>■ x. Figure 1 illustrates what 

happens close to t — 0. 

(b) In a more general case, ie for arbitrary values of our parameters, the configuration 
represents a multi-soliton solution that is localized close to the centre of the mass and 
deforms as the time passes. As an example, let us present a typical case by taking f(z) = z, 
g(z) = 2z, h(z) = z 2 and h\{z) = 3. In this time-dependent solution, which corresponds 
to a configuration of three solitons, the energy density of the field, for large (negative) t, 
is in the shape of a ring which, as time passes, deforms to three peaks and then expands 
again to a ring. Figure 2 presents pictures of the corresponding energy densities at some 
representative values of time. 

Returning to the scattering of our soliton-like structures we have seen that an embed- 
ding of the £77(2) model provides an example of such a scattering. A further example is 



provided by a genuine SU (3) case when we take 

f = p 2 (z), g = V2p(z), (28) 
where p(z) is an arbitrary function. 

One example is presented in Figure 3, where we have chosen p = z, h = z 3 and 
hi = 0. This configuration can be seen to consist of two static solitons at the centre- 
of-mass accompanied by two solitons accelerating towards the centre, scattering at right 
angles, and then decelerating as they separate. 

V. CONCLUSIONS AND FURTHER COMMENTS 

In this paper we have shown how to adapt Uhlenbeck's construction of uniton solutions 
of the chiral model in (2+0) dimensions to the, time-dependent, uniton-like solutions of 
the modified chiral model in (2+1) dimensions. It seems likely that there are many more 
interesting solutions still to be found. One could, for example, investigate the existence 
of time-dependent solutions based on other static ones constructed by Uhlenbeck. 

As for static fields both models are the same, all Uhlenbeck's solutions are also solutions 
of the modified chiral model. However, at the level of two unitons and beyond, our 
construction has given us genuine time-dependent solutions of the modified chiral model. 
Looking at some examples we have seen that, in some cases, their field configurations 
can be thought of as representing scatterings of some soliton-like objects, sometimes 
accompanied also by static structures. The total energy appears to be quantized - thus 
showing that during their scattering the soliton-like structures must change their size. 

We have not succeeded yet in gaining the full understanding of various properties of 
our solutions. These problems are currently being investigated. 
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Figure Captions. 

Figure 1: Energy density of the embedding of SU(2) solitons, at increasing times. 

Figure 2: Energy density at increasing times, for a three ring-shaped solitons. 

Figure 3: Energy density at increasing times showing a 90° scattering between pure 
SU(3) solitons. 
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